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Active particles contain internal degrees of freedom with the ability to take in and dissipate energy 
and, in the process, execute systematic movement. Examples include all living organisms and their 
motile constituents such as molecular motors. This article reviews recent progress in applying the 
principles of nonequilibrium statistical mechanics and hydrodynamics to form a systematic theory 
of the behaviour of collections of active particles - active matter - with only minimal regard to 
microscopic details. A unified view of the many kinds of active matter is presented, encompassing 
not only living systems but inanimate analogues. Theory and experiment are discussed side by 
side. This review is to appear in volume 1 of the Annual Review of Condensed Matter 
Physics in July 2010 and is posted here with permission from that journal 



I. INTRODUCTION 



Systems of Interest: Flocks, Rods, and the 
Cytoskeleton 



Physics enters biology in two broad and overlapping 
areas - information and mechanics [l|, 0. This re- 
view is about the mechanics, statistical and otherwise, 
of living matter. This is not a review of soft-matter 
physics in a biological setting. The aspect of biological 
matter of interest here is the ability to transduce free 
energy into systematic movement. This property is 
the defining characteristic of active matter, and my in- 
terest is in the unique mechanical properties that arise 
as a consequence of such active processes [1] . This re- 
view focusses on the collective behaviour of systems 
with a large number of constituents, and will thus use 
the ideas of condensed matter and statistical physics. 
I will not discuss the physics of molecular motors or 
the related and vast exclusion-process family of prob- 
lems [J, [6|? ] . The aim is to complement and update 
the perspectives of earlier reviews @, [1] and to present 
the achievements and limitations of this rapidly ad- 
vancing subfield. I begin with the realizations of ac- 
tive matter to be covered. 



'Electronic address: sriram@physics.iisc.cmct. in 
tAlso with CMTU, JNCASR, Bangalore 560 064 



Flocks: theoretical, real, and imitation 



It is natural for a condensed matter physicist to re- 
gard a coherently moving flock of birds, beasts or 
bacteria as an orientationally ordered phase of liv- 
ing matter. This idea was first implemented 0, H3] 
in computer models of interacting particles moving 
with a fixed speed and trying to align their velocity 
vectors parallel to those of their neighbours, in the 
presence of noise. To a student of magnetism these 
are itinerant, classical, ferromagnetically interacting 
continuous spins that move in the direction in which 
they point, which is what makes flocks different from 
magnets. The models showed a nonequilibrium phase 
transition from a disordered state to a flock with long- 
range order [Ilil - fl3j in the particle velocities as the 
noise strength was decreased or the concentration of 
particles raised. The nature of order and fluctuations 
in a flock, the character of the transition to a flock, 
and the flocking of active particles with an axis but no 
polarity are among the issues this review deals with. 

Amazingly, the physics of flocking can be imitated 
by a collection of rods lying on a horizontal surface, 
agitated vertically [l4j]. Indeed, vibrated monolay- 
ers of macroscopic grains have provided some of the 
most fruitful realizations of active matter. This review 
will explain briefly the physics behind underlying this 
analogy and highlight experiments that exploit it (f5j 



Flocking in a fluid medium introduces physics absent 
in the simplest flocking models: each swimming crea- 



ture generates fluid flow which moves and reorients 
other swimmers far away. Including this interaction 
leads [l6| to a modified liquid-crystal hydrodynamics 
[l7^ . in which the constituent particles carry perma- 
nent stresses that stir the fluid. The implications of 
this interplay of self-propelling activity and fluid flow 
for order and macroscopic rheology are an important 
part of this review. 



2. The cytoskeleton as an active gel 



The physics of the flocking of organisms in a fluid 
reappears at a siz&cellular scale in the cytoskeleton 
@) Ell > the polymeric scaffolding that governs trans- 
port, adhesion, movement and division in the liv- 
ing cell. Two nonequilibrium processes drive the cy- 
toskeleton: the ATP-assisted polymerization and de- 
polymerization known as treadmilling, which I will not 
discuss; and contractility, illustrated in Fig. [TJ which 
arises from the ATP-driven movement of motor pro- 
teins in specific directions on biofilamcnts, and is cen- 
tral to this review. The cytoskeleton is thus a suspen- 
sion of filaments endowed with active internal forces. 
Moreover, there are natural mechanisms that promote 
the alignment of neighbouring filaments, through ex- 
cluded volume as well as activity. We must there- 
fore allow for the possibility of orientational order. It 
is then not surprising that the hydrodynamic equa- 
tions obtained in the active-gel description (l9l. l2dj] 
of the cytoskeleton have precisely the same form as 
those [16| for collections of swimming organisms, ig- 
noring complications such as permanent crosslinking. 
A condensed-matter physicist accustomed to symme- 
try arguments should find it reassuring that identi- 
cal hydrodynamic descriptions apply for a 10 km fish 
shoal and a 10 /j,m cytoskeleton. 



B. Viewpoint: active systems as material 



The viewpoint of this review is that living matter can 
fruitfully be regarded as a kind of material J21j and 
studied using the tools of condensed matter physics 
and statistical mechanics; that there is a practical way 
to encode into such a description those features of the 
living state that are relevant to materials science; and 
that the results of such an endeavour will help us bet- 
ter understand, control and perhaps mimic active cel- 
lular matter. A community of respectable size has 
grown around this activity, approaching the problem 
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FIG. 1: A cluster of motors with heads on both sides, 
exerting contractile forces on the cytoskeletal network. 

at different length scales and with a variety of tech- 
niques. What do they, and we, hope to gain from this 
enterprise? First, active matter is condensed matter in 
a fundamentally new nonequilibrium regime: (i) The 
energy input takes place directly at the scale of each 
active particle, and is thus homogeneously distributed 
through the bulk of the system, unlike sheared fluids 
or three-dimensional bulk granular matter, where the 
forcing is applied at the boundaries, (ii) Self-propelled 
motion, unlike sedimentation, is force-free: the forces 
that particle and fluid exert on each other cancel, (iii) 
The direction of self-propelled motion is set by the 
orientation of the particle itself, not fixed by an exter- 
nal field. Indeed, these can be taken as the practical 
defining properties of active matter. A comprehensive 
theory of this ubiquitous type of condensed matter 
is a natural imperative for the physicist, and should 
yield a catalogue of the generic behaviours, such as 
nonequilibrium phases and phase transitions, the na- 
ture of correlations and response, and characteristic 
instabilities. Second, therefore, the generic tenden- 
cies emerging from the theory of active matter, unless 
suppressed by specific mechanisms, must arise in vivo, 
which is why biologists should care about it. Last, 
if we can understand active matter, perhaps we can 
manufacture faithful imitations of it in chemomechan- 
ical systems without components of biological origin 

The reader should keep in mind that theories of active 
matter were formulated not in response to a specific 
puzzle posed by experiments but rather to incorpo- 
rate living, metabolizing, spontaneously moving mat- 
ter into the condensed-matter fold. This was done 
through minimal models whose consequences are rel- 
atively easy for theoreticians to work out. Natural 



realizations of living matter are far from minimal, so 
comparisons of active-matter theory with experiment 
are likely to be qualitative until well-controlled model 
systems are devised. 



C. Prehistory 



Although this is not a historical review, it must be 
stated here that the active-matter idea is not new. In 
a prescient article, Finlayson and Scriven [24j argued 
that biological matter could display hydrodynamic in- 
stabilities driven by stresses arising from metabolic ac- 
tivity in the bulk of a fluid. Unlike in this review, they 
build active stress-tensor contributions from gradients 
of existing scalars such as concentration and temper- 
ature. They point out that the field poses "challenges 
to theoretical rheology" and is "virgin territory for 
experimental research in the physics of complex flu- 
ids". They conclude that "continuum., analysis... of 
instability by active stress is probably .... destined 
to play a large part in understanding some of the en- 
gines operating at the cellular level in living systems" . 
The developments surveyed in this article prove them 
right. 



D. Structure of this article 



The body of this review consists of several sections 
each treating a specific example of an active-matter 
system. In Section |TT] the simplest models of moving 
flocks, without solvent flow, are introduced and some 
key properties discussed. Flocks without macroscopic 
polarity and hence no migration velocity are surveyed 
in section irm The theoretical framework for collective 
self-propulsion in a fluid medium is found in Section 
IIV1 with a discussion of general symmetry-based ap- 
proaches, and the active mechanics of the cytoskeleton 
in Section IIVF1 Non-biological self-driven particles 
are the subject of section [Vj The article closes with 
Section I VII which summarizes the achievements and 
limitations of the approach adopted here to the mod- 
elling of active matter and attempts to set a course 
for future experiments and theory in this vital and 
rapidly evolving field. 



II. POLAR FLOCKS ON A SUBSTRATE 



Arrows are polar objects, uniform cylinders are apo- 
lar. A flock moving in one direction has polar order, 
whereas equal numbers of ants moving from south to 
north and from north to south are macrscopically ap- 
olar even though the individual ants are polar. We 
begin by discussing polar ordered states. 



A. The Vicsek model 



In the simplest models [9|,|l0( of ordering transitions of 
active particles, the medium through which the parti- 
cles move is treated as an inert substrate. Each parti- 
cle carries a velocity vector of fixed magnitude, which 
it rotates to point parallel to the mean of its neigh- 
bours' velocities, with an angular tolerance rj. It then 
takes a small step in the direction of the updated vec- 
tor. Like the continuous-spin magnets that they re- 
semble, the Vicsek family of models display a well- 
defined phase transition from a disordered phase to a 
coherent flock as rj is decreased or the number de nsity 
is increased [To|, [H, [H| . Unlike planar spins [27l.l28j. 
flocking models show true long-range order even in 
two dimensions. How they do so was shown by Toner 
and Tu @, El, [Hj], whose continuum field-theoretic 
approach, which amounts to a coarse-graining of the 
Vicsek rules, we sketch below. 



B. Toner- Tu field theory 



A coarse-grained dynamical description works with 
"slow variables" [29[ , whose relaxation times increase 
unboundedly with increasing wavelength. For the Vic- 
sek model these are [III E3 the number density field c, 
since total number is conserved, and the velocity field 
p, whose fluctuations transverse to the mean order- 
ing direction are the spin- waves or Nambu-Goldstone 
(30L l3l| modes of the flock f32|. With the slow vari- 
ables in hand, Toner and Tu [ill. [l2| write down the 
equations of motion, incorporating the physics of the 
Vicsek model, on general grounds of symmetry. 

First, the dynamics of p is described in the equation 
below: 

<9 t p + ApVp + ... = (a-/3p-p)p + rVVp-VP(c) + f 

(1) 



Ignoring all gradients in Equation [TJ a phase transi- 
tion from the is otrop ic state p = to an ordered flock 
with | (p) | ~ y/a/P, spontaneously breaking rotation 
invariance, should occur as the control parameter a is 
tuned from negative to sufficiently large positive val- 
ues; T, in general a tensor, controls the elastic [33| 
restoring forces to distortions in the ordered phase. 
The random error 77 of the Vicsek update rule is sum- 
marized in a nonconserving, Gaussian noise f , uncor- 
rected in space and time. The crucial difference from 
traditional models of dynamic critical phenomena 34] 
lies in the Navier-Stokes-like A term [35|], which says 
that distortions in p are advected by p, because p is 
not only an order parameter but a velocity. P(c), a 
general increasing function of concentration, then em- 
bodies the equation of state for the pressure. Second, 
c evolves through the continuity equation 

d t c + V-cp = 0, (2) 

because p is the velocity field of the particles. Equa- 
tions [TJ and [5] can be viewed as the dynamics of a 
fluid with a preferred speed relative to a background 
medium, or of a magnet whose spin is a velocity. In the 
second interpretation, the nonequilibrium character of 
the Toner- Tu model enters only through the advective 
A term in Equation QJ and the current in Equation [2j 
The right-hand side of Equation [TJ can be written as 
—SF/Sp plus noise, with a free-energy functional 

F[p] = J d d r[-|p.p + ^(p.p) 2 + ^VpVp+p-VP(c)]| 

(3) 

of Ginzburg-Landau form, in which the "pressure gra- 
dient" VP(c) appears as an orienting field for p. In 
the absence of A, the average dynamics is downhill in 
F. A flock governed by Equations [TJO is a strange 
blend of magnet and fluid. Let us rediscover some of 
its unique properties [ll|, UH , all consequences of the 
fact that the order parameter is itself a velocity. 



C. Sound modes, giant fluctuations, and 2d 
long-range order 



Sound modes: Let us linearize Equations [TJ [2] about a 
uniform ordered state with |(p)| = po, c = cq, denote 
directions along and transverse to (p) by || and _L, 
and see how small disturbances 5c, Sp = (5p±,5p\\) 
travel. As in Heisenberg or XY magnets [33], we can 
see from Equation [TJ that <5p|| relaxes rapidly, so that 
on long timescales Sp ~ 5p±. If we take A — 1 in 
Equation [TJ a shift to a frame moving with speed po 
eliminates advective terms from Equations [TJ and [2] 
in this linearized analysis without losing the essential 



physics. To leading order in gradients, and ignoring 
the noise, Equations [TJ and d] then become 

d t 5c= -c V ± -8p ± , (4) 

dt6p± = -P'(c )V±5c (5) 

leading to propagating modes with frequency 

u = ±^/c P'{c )q ± (6) 

at wavevector q = (q±,qn). No, we haven't rediscov- 
ered normal sound waves. Long- wavelength sound in a 
fluid propagates because of momentum conservation, 
which a flock lacks. The "sound modes" in a flock are 
a consequence of spontaneously broken rotation in- 
variance: they propagate at long wavelength because 
5p± is not damped at zeroth order in wavenumber. 

Giant number fluctuations: Sp± is the Nambu- 
Goldstone mode [3(3, HU of broken rotation invariance, 
easily excited and slow to decay at small wavenumber 
q. Hence when noise is included the steady-state vari- 
ance (|5p± q | 2 ) of the qth fourier component should 
diverge at small wavenumber q, like spin waves in 
an XY magnet [37j or director fluctuations in a ne- 
matic [13] • But Equations [4][5] imply that, for fre- 
quencies and wavenumbers related by Equation (BJ 
8c ~ y / Co/P'(co)8p±. Therefore, the variance of 5c 
should also diverge at small q. In real space this 
means that for regions with N particles on average, 
the variance in the number grows faster than N: a 
flock should display giant number fluctuations. 

Long-range order in 2d: Equations [JJ and [5] when lin- 
earized yield (|6pj_ q | 2 ) ~ l/q 2 - Taken literally this 
would mean only quasi-long-range order in dimension 
d = 2, by analogy with the XY model. However, non- 
linearities in Equations [TJ and [4] are stongly relevant 
in d = 2. and lead to (|<5p_L q | 2 ) diverging more slowly 
than ~ 1/g 2 for most directions of q, thus preserving 
long-range order [III Qjj]. Here is a qualitative ex- 
planation of how this happens. Consider a flock with 
(p) = Po% at a given instant. Because p is a velocity, a 
long wavelength fluctuation 5p± will allow creatures 
in the flock to exert their orienting influence on re- 
gions out to a distance 5p±t in a time t. This effective 
enhancement of the range of influence ultimately leads 
[TTl - fT3| to a suppression of orientational fluctuations 
on large scales, and thus to long-range order in d = 2. 
The effect is absent in d = 1 , where rotations are a dis- 
crete symmetry and there is no 5p±. We will return 
to this point in section Hi D 31 

I have not discussed microscopic derivations of the 
Toner- Tu equations f38l l39ll. dynamics near the on- 
set of polar ordering [40, [4J| , and the application of 



flocking models [42], |43[ to motor-microtubule extracts 
0- 



D. Tests of vector flocking models 



Two types of tests of the physicist's approach to flock- 
ing are required. Experiments on real flocks are of 
course the real test. However, it is equally impor- 
tant to check the predictions of coarse-grained theories 
against computer experiments on microscopic models. 
We discuss both below, although measurements of the 
phase-transition variety on real flocks are scarce. 



transition in the limit of infinite system size, as occur§ 
in Vicsek models in dimension d > 2. Recall that the 
Id Ising model fails to order because of the prolifer- 
ation of kinks 37 1. Does some magic suppress kinks, 



i.e., velocity flips, in Id flocks? Raymond and Evans 
[Hol ] estimate the lifetime r(L) of a Id flock of length L 
as the mean time for the appearance of a kink. They 
show that if t(L) is to grow at least as a power of 
L, the number of particles per site must grow, ar- 
tificially, at least as logL. Accepting the constraint 
of a finite number density yields a finite lifetime and 
correlation length for the would-be flock. Ultimately, 
nothing saves Id flocks from the fate of the Id Ising 
model because flocking in Id breaks a discrete symme- 
try and is thus not accompanied by Nambu-Goldstone 
modes whose long range is responsible for the rescue 
act in d = 2. 



1. Scaling in the ordered phase - simulations 



Quantitative agreement has been found between nu- 
merical "experiments" on microscopic models and the 
predictions of the coarse-grained theory, including 
long-range order in d = 2, the form of the propagating 
modes, anomalous density fluctuations and superdif- 
fusion of tagged particles 0, El, [H[ . I am not aware 
of any laboratory or field-based attempts to measure 
ordered-phase correlation functions for flocks on a sub- 
strate. 



2. Nature of the flocking transition ~ simulations 



Although some studies [T3, EH, IH| claim to see a con- 
tinuous onset of the ordered phase, as a mean-field 
treatment of EquationQ]or analogy to continuous-spin 
magnets would suggest, very large-scale numerical 
studies of the Vicsek model and its relatives show that 
the transition is in general a discontinuous one, char- 
acterized by a complicated coexistence [H, US E3j- 
See also Ref. gg. 



3. On one- dimensional flocks 



In numerical studies of the Vicsek model in one space 
dimension, long domains of coherently moving parti- 
cles do appear at low noise and high density. It is 
occasionally claimed IH| that this is a true phase 



Nevertheless, lowering the noise strength or raising 
the density at fixed L will eventually yield a system 
with correlation length larger than L, and a bistable 
time-series of the flocking order parameter, i.e., a finite 
transient flock with a large lifetime. Precisely this 
phenomenon has been seen in experiments, to which 
we turn next. 



4- Real experiments on flocks on a substrate 



1. Couzin 51] studies the one-dimensional flock- 
ing of desert locusts on an annular track. The 
onset of coherent motion seen in the experi- 
ments as the number of locusts is increased, and 
the bistable time-series switching between clock- 
wise and counterclockwise circuits of the track 
are well described by the Id Vicsek model [Toj . 
to which the paper makes a comparison. One 
must not of course conclude - see Reference 
and section III D 31 - that there is a true phase 
transition here. All the same, the experiment 
demonstrates convincingly the relevance of sim- 
ple flocking models to an understanding the be- 
haviour of real organisms. 

2. Szabo et al. 52] created a polar flock in a petri 
dish by using motile keratocytes extracted from 
fish scales. As the cells divide and the areal 
density of cells crosses a threshold, there is a 
well-defined onset of a state of macroscopic or- 
der in the orientations and velocities of the ker- 
atocytes, that can plausibly be interpreted as a 
flocking phase transition, which should be gov- 
erned by the local physics of the models we have 
been discussing, as the cells interact only by con- 



Equal friction: High friction Low friction 

symmetric propulsion 

FIG. 2: A rod landing after being tossed up will in general 
be impelled away from the end that makes first contact 
with the surface. If the two ends differ in weight, geometry 
or friction, the rod will be propelled towards one end. 

tact. The data quality makes it impossible to 
draw conclusions regarding the order of the tran- 
sition. The authors make further inferences on 
the transition from a computational model, not 
from data. Despite limitations the experiment is 
valuable as a realization of the flocking transi- 
tion, with potential for higher-precision studies 
including finite-size scaling. 

3. There is a surprising link between vibrat ed g ran- 
ular monolayers and self-propulsion 0, 02 HH, 
HI, HH ■ Polar rods [23[ , placed on a vertically 
vibrated horizontal linear track, move on aver- 
age with small end forward. Systematic uni- 
directional motion of apolar rods [53J was ef- 
fected by maintaining them permanently tilted, 
by overfilling an annular channel with rods. The 
propulsion mechanism [53 - i56j is illustrated in 
the cartoon Fig. [2j The reader is encouraged to 
check that this system satisfies all three defining 
properties of active matter in section II Bl Pub- 
lished studies so far [57} of polar rods at high 
coverage on a vibrated surface report order and 
anomalous density fluctuations on a local scale, 
and an abrupt flight to the boundary when the 
vibration amplitude crosses a threshold, but not 
a macroscopically ordered flock. The case of ap- 
olar rods (la, [58| is discussed in section IIIIB 21 



III. FLOCKS THAT GO NOWHERE ON A 
SUBSTRATE 



A. The origin of giant number fluctuations in an 
active nematic 



At least one example is known [59] of living, pulsat- 
ing elongated cells forming a state with spontaneous, 
uniaxial orientational order with a macroscopic axis n 
with n > — n symmetry. This is the spatial symmetry 
of nematic liquid crystals, and it guarantees that the 
mean macroscopic velocity of such a "flock" is zero. 
And yet active nematics are profoundly different from 




ASYMMETRY LEADS TO CURRENT 
J x cx d z 6 J z <x d x 8 

FIG. 3: In an apolar nematic, curvature confers local po- 
larity and hence local motion. 



dead ones. In particular, here as in polar flocks, giant 
number fluctuations are predicted [l4j - even though 
the order parameter for the nematic is not a velocity. 
Below is an explanation of how it works. 

In the presence of noise, long wavelength fluctuations 
of the axis of orientation - the broken-symmetry mode 
- are abundant. Working in two dimensions for sim- 
plicity, where distortions are described by a single an- 
gle 9, we see (Fig. [3]) that bend or splay produces a 
polar configuration. The absence of time-reversal in- 
variance in a driven state, together with the polarity, 
means [6(| that particles in the curved region must 
drift in a direction defined by the curvature as pic- 
tured in Fig. |31 leading to a current of particles J 
with J x cx d z 9, J z oc d x 9 in two dimensions. In steady 
state, this active flux must be balanced on average by 
restoring diffusive currents ~ Vc. Thus fluctuations 
in c and 9 must be of the same order. Orientational 
fluctuations, being a broken-symmetry mode, should 
have a variance of order 1/q 2 at wavenumber q. There- 
fore, number density fluctuations (|c q | 2 ) ~ 1/q 2 as 
well which means the number fluctuations in regions 
containing N particles on average have a standard de- 
viation 



AN cx N 1/2+1,d ~ N for d = 2 (7) 



for large N, Unlike for polar flocks [TT|, these re- 
sults of the linearized theory for apolar systems sur- 
vive quantitatively when nonlinearities are included 

Mm- 



B. Tests of apolar flock theories: experiments 
and simulations 



1. Particle models for active apolar order 



The first test of the existence of giant number fluc- 
tuations (GNFs) in active nematics was a computer 
experiment on a clever apolar generalization [62J of 
the two-dimensional Vicsek model described in section 
III A[ Each particle is assigned an axis, rather than a 
unit vector, which it aligns parallel to the mean of 
its neighbours' axes, subject to a small angular noise. 
The particle then takes a small step preferentially 
along the axis, forward or back. Since the orienta- 
tion variable and the moves are all defined modulo 7r, 
the model has no polarity. As noise or mean inter- 
particle separation is decreased, a continuous transi- 
tion is seen from an isotropic state to a phase with 
quasilong-range nematic order of the particle axes as 
in equilibrium fTA l63j. However, the density fluctua- 
tions are giant: AA~ grew far more rapidly than N and 
could be fitted to a form ay/N + bN in precise agree- 
ment with the predictions [14| whose derivation was 
sketched above. Real-space snapshots showed strong 
transient inhomogeneity, with a banded structure. If 
the anisotropic, detailed-balance-violating move is re- 
placed by an isotropically distributed random step, 
the nematic phase survives, but with normal number 
fluctuations, because detailed balance now holds. S. 
Mishra [39| has shown that coarse-graining the model 
of Chate et al. yields the equations of motion 

of Ref. Q, with the active current as in section 
MI Al Further numerical studies on a particle model 
64] show a close connection between the giant number 
fluctuations of an active nematic and the phe nomenon 
of fluctuation-dominated phase ordering 65] . 



2. Experiments on apolar flocks on a substrate 



1. Melanocyte nematics: Living melanocytes, the 
cells that spread pigment in skin, have been 
shown [59l | to form apolar, nematic order in vitro 
as concentration is increased. The resulting 
2d nematic phase has large regions of oriented 
cells with occasional point topological defects of 
strength —1/2 which reinforce the nematic inter- 
pretation [5^] . Although the authors of [E^] 
emphasize that their systems are very far from 
thermal equilibrium the analysis they bring to 
bear on the problem, especially in the apolar 
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FIG. 4: Snapshot from the active nematic phase of 
vibrated copper-wire segments. Figure credit: Vijay 
Narayan 

case, is largely inherited from the equilibrium 
theory of liquid crystals. The focus in 59] is 
on modelling cell orientations, in a population- 
averaged description which cannot resolve spa- 
tial variations in the density. Consistent with 
this limitation, their experiments do not look at 
the behaviour of density fluctuations. The pur- 
pose of this discussion is to urge experimenters 
to take another look at the melanocyte system, 
especially as the strange properties of active ne- 
matics m m n remain to be checked in a 
biological experiment. 

2. Active apolar states in granular matter: 

Narayan et al. [HI, [H[ carried out a systematic 
study of liquid-crystalline order in a vertically 
agitated monolayer of copper-wire segments. In 
particular, they measured number density fluc- 
tuations in the nematic and the isotropic phases. 
Their findings fl5l] confirmed the predictions 
[Tij of giant number fluctuations in the nematic 
phase. A typical nematic configuration and a 
plot of standard deviation vs mean number are 
shown in Figs. [4] and [5] An additional con- 
sequence of the theory that the autocorrelation 
of the local density should decay as — log t over 
a large time range [l4| was also confirmed (l5j . 
Lastly, single-particle tracking showed that the 
rods moved preferentially along their length just 
like the apolar flockers in the model of Chate et 
al. 

Other experiments of interest, on chiral struc- 
tures in active granular matter, include the work 
of Tsai et al. [66| and of Blair et al. [53| . 
The measurements on granular matter reported 
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FIG. 5: Giant number fluctuations in experiments on the 
copper-wire active nematic: standard deviation AN scaled 
by \HV grows with mean number N. Inset: no such growth 
in the isotropic phase. Figure credit: Vijay Narayan 



above are among the few quantitative tests of 
active-matter theories. That they are not car- 
ried out on living systems only serves to rein- 
force the universality of the approach reviewed 
in these pages. 



IV. SWIMMING ORGANISMS AND THE 
CYTOSKELETON 



A. The equations of active hydrodynamics 



The collective motion of swimming organisms can be 
understood by constructing the coupled dynamical 
equations for the swimmer concentrations and orien- 
tations, and a generalized Navier- Stokes equation for 
the velocity field u of the suspension. Self-propelling 
activity enters via force densities in the Navier-Stokes 
equation. These have no monopole moment, since the 
mutual forces of swimmer and fluid cancel by New- 
ton's 3rd Law (67j . The minimal model of an active 
particle in a fluid is therefore a permanent force dipole 
(see Fig. [6]), whose strength W is positive for exten- 
sile and negative for contractile swimmers [See also 
the modelling of ion pumps in a model membrane by 
Manneville et al. j68j.] Associating a local orienta- 
tional order parameter field p(r, i) with the vectorial 
orientation of the swimmers in a coarse-graining cell 
around r at time t, a straightforward calculation 
shows, to leading order in a gradient expasion, that 




FIG. 6: The flows around a swimming bacterium (top) 
and an algal cell (bottom): bacteria are extensile, and 
algae contractile, force dipoles. 

a collection of active particles at concentration c(r, t) 
has an active force density F a (r, t) = — V • a a with an 
active stress 

a a = Wc(r,*)pp 

W „ 

= Wc(r,t)Q(r,t) + —cp 2 \, (8) 

where I is the unit tensor, and Q, the trace-free part of 
pp, determines the deviatoric stress which is all that 
concerns us in the incompressible limit V • u = 0. Ac- 
tivity thus produces flow, and flow in turn reorients 
the principal axes of rodlike particles preferentially 
along the extensional axis (and platelike along the 
compressional axis), as it would in non-living liquid 
crystals as well. The resulting equations of motion for 
p are obtained by adding flow-orientation couplings 
to Equation [TJ 

SF 

Ap + Ap- Vp + .... = 7 A p- — +f, (9) 

dp 

where D t is the time-derivative in a frame locally co- 
moving and corotating with the fluid, the coefficient 
7 governs the orienting effect of the extensional part 
of the flow through the term in A = [Vu + (Vu) T ]/2, 
and the remaining terms are as defined for Equation 
[Tl F being the free-energy of Equation [3] Inertia, vis- 
cosity 77, elastic forces from Equation [3l and activity 
compete to determine the dynamics of the suspension 
velocity field: 

p(<9 t +u-V)u = -V-(cr a +cr' l )-ryV 2 u-Vn; Vu = 0. 

(10) 

with incompressibility V • u = enforced by the pres- 
sure II. In Equation [TU1 a h is the contribution to the 




FIG. 7: The heart of the generic instability of ordered 
active filaments: long- wavelength splay produces shear 
flows that further distort a row of parallel contractile force 
dipoles (left); similarly bend disrupts extensile filaments 
(right). 



stress from Equation [3] for the free-energy functional 
[l7j . Lastly, the current in the continuity equation 
for the concentration, in the lab frame, is (p + u)c, 
because we can simply identify the orientation field p 
with the active-particle velocity with respect to the 
fluid. 

We turn next to the dramatic consequences of this 
dynamics. 



B. The instability of orientationally ordered 
active suspensions 



tainer with smallest dimension L, active stresses corrr^ 
pete (70l]with orientational elasticity from Equation^ 
resulting in a threshold Wcq ~ K/L 2 for the instabil- 
ity. Far beyond threshold, scaling at wavenumber q 
yields a growth rate (qL) 2 Wco/r] for qL <C 1, crossing 
over to Wco/r] for qL ^> 1. The detailed form of the 
instability depends on the experimental geometry. In 
Ref. 69], for example, it appears as an active version 
of the Freedericksz transition of nematic liquid crys- 
tals [13] ■ Novel instability mechanisms recently noted 
in polar systems include oscillations and banding if a 
concentration field is included [72[ , and travelling un- 
dulations if a deformable free surface is present f73j . 
Aranson et al. [74[ study free-standing active thin 
films and introduce ad hoc a propulsive force simply 
proportional to the thickness-averaged p. Such a term 
can be obtained from the treatment of this section and 
Ref. if p is neither parallel nor perpendicular to 
the free surface. The stabilizin g eff ect of shear has 
been discussed by Muhuri et al. |75j |. 

There is detailed evidence for the basic instability in 
numerical studies of particles self-propelled through a 
fluid (76Tj83j , who also point out a variety of instabil- 
ities in the isotropic phase. Studies of the PDEs of 
active hydrodynamics include Refs. (85f ? - 489j who 
show that the linear instability leads to complex flow 
patterns and turbulence driven not by Reynolds num- 
ber but an active Ericksen number (Wco/r])T where r 
is the collective relaxation time of the Q tensor. 



Active suspensions with uniform orientational or- 
der are in a state of permanent uniaxial tension or 
compression and are therefore intrinsically unstable 
[H, |6^, : An individual contractile particle (Fig. [5]) 
pulls fluid in from both ends along its main axis. In 
a perfectly ordered, unbounded, parallel collection of 
contractile particles, these self-generated flows cancel. 
A long- wavelength splay fluctuation disrupts this deli- 
cate balance, resulting in flows upward on one side and 
downward on the other, as shown in Fig. [7] The ensu- 
ing shear will amplify the rightward tilt of the middle 
portion, hence the instability. The same argument, 
mutatis mutandis, implies a bending instability of ex- 
tensile filaments [7l| . A detailed solution [Tgj, [6^, [7(| 
of Equations M ITU1 bears out this pictorial argument. 

Note that we've assumed instant, global response of 
the fluid flow, which amounts to the steady Stoke- 
sian approximation in which inertia and acceleration 
are ignored in Equation 1101 In that limit, in an 
unbounded system with mean concentration Co, the 
growth rate of the instability must be proportional 
to the only available inverse timescale, namely, the 
ratio Wco/rj of active stress to viscosity. In a con- 



For small but nonzero Reynolds number Re on the 
length scale a of a particle, the instability arises for 
wavenumber q with 

Re a < qa -C 1, (11) 

a range that clearly does not exist for swimmers of 
macroscopic size and speed. 



C. The effect of activity on viscosity 



The enhancement or reduction of viscosity by activity 
[90( | in a suspension in the isotropic phase is one of 
the most robust predictions of active hydrodynamics. 
The mechanism is readily understood, pictorially (see 
Fig. or verbally: consider a collection of filaments 
endowed with contractile or extensile force dipoles, in 
the isotropic phase but at high enough concentration 
that collective orientational relaxation is slow. A mod- 
est imposed shear flow will then produce appreciable 
alignment. Contractile filaments will then pull back 
on, and extensile filaments push out on, the flow that 





FIG. 8: How activity modifies viscosity: Shear orients fil- 
aments; the permanent force dipoles pull back, if contrac- 
tile, and push out, if extensile, on the flow. 



oriented them. The result is a higher stress/rate ratio, 
that is, a viscosity, for contractile filaments, and lower 
for extensile filaments, than would be found for dead 
particles of the same shape and at the same concentra- 
tion. Upon tuning a parameter such as concentration 
to increase the collective orientational relaxation time 
r, a suspension of contractile filaments should show 
strong viscoelasticity, like an equilibrium system ap- 
proaching a glass transition [9fJ, l£l| and extensile fila- 
ments should show enhanced unstable shear-thinning. 
Numerical studies [§2[ find a complex active rheology 
including the predicted [90, l9l| glassy enhancement of 
viscosities for contractile systems, and shed valuable 
light on the role of boundary conditions. Microscopic 
models for filaments and motors 1911 and swimmers 
[931 ) recover the predictions of Ref . [90| ■ 



Active noise in the ordered phase is predicted jl04j 18 
lead to diffusivities that depend on sample thickness, 
through a mechanism related to the giant fluctuations 
[Pll - fTl llt| discussed in sections HT1 and MI Al 



E. Fish schools 



Schools of fish [1051 ] appear to be sponta neous ly or- 
dered phases at high Reynolds number [l06| and, 
from Equation 111! evade the instability of the ori- 
ented state of low-Re self-propelled particles. As a 
first step towards applying active hydrodynamics to 
macroscopic swimmers, Ref. [lj| studied linear per- 
turbations about an ordered steady state of uniform 
number density and propulsion velocity relative to a 
quiescent fluid. Ignoring inertia and retaining active 
stresses and the Navier- Stokes acceleration term in 
Equation llOl they found that all modes were propaga- 
tive, with highly anisotropic speeds of order ^Wco/p, 
and showed that the inclusion of fluid flow did not 
eliminate the giant number fluctuations in polar flocks 
predicted by Toner and Tu fllj ] . It would seem worth- 
while to test these predictions on fish schools. See Ref. 
107], to which we return in section HV G 31 below. . 



D. Fluctuating active hydrodynamics 



The enhancement of biological noise by active pro- 
cesses was studied some years ago in the context of 
membranes with pumps [68|, 1941-98] . A natural expla- 
nation for the enhanced noise temperature of swim- 
ming bacteria [§9[ was found in the active hydrody- 
namic framework by including random forces, torques 
and currents (l6l . |9G ] which are consistent with con- 
servation laws but unrelated a priori to transport co- 
efficients, because these are syste ms far fro m thermal 
equilibrium. Lau and colleagues [lOOl llOll ] model the 
microrheology of bacterial suspensions in the active 
hydrodynamic framework with noise. They relate the 
collective noise to the tumbling of bacteria, and show 
that Equation [3] implies a t~ x / 2 decay of the stress 
autocorrelation for times shorter than the collective 
orientational relaxa tion time r. They also show, as 
do Underhill et al. |102| |. that equal- time correlations 
of the suspension velocity field on intermediate length 
scales r should decay as 1 jr. R elated predictions arise 
in work by Golestanian jl03| in the a priori differ- 
ent context of catalytically self-propelled colloids [22[ . 



F. Active hydrodynamics of the cytoskeleton 



The cytoskeleton with its motors, on timescales long 
compared to the release time of transient crosslinks, 
is also an active fluid susp ension of polar orientable 
objects [1, d HO, Efil 109] with macroscopically con- 
tractile stresses along the filament axes [l 10j . 



1. The active-gel picture 



Hydrodynamic equations for active gels [8l. [T9l[20lllO9| 
were formulated with a view to modelling the cy- 
toskeleton, mainly on long timescales when it is fluid. 
The description turns out to be precisely the same as 
that of Ref. [16], but the const ruct ion uses linear ir- 
reversible thermodynamics p9l . Illl| for a nonequilib- 
rium steady state near equilibrium. A linear relation 
is assumed [h| between "fluxes" (including the stress 
o~ij and the rate r of ATP hydrolysis) and "forces" 
including the velocity gradient V^ity and the chemical 



potential difference A/i between ATP and its reac- 
tion products, allowing for a nondiagonal matrix of 
kinetic coefficients. The existence of the polar order 
parameter p makes it possible to construct a scalar 
CPiPj^iUj contribution to r, linear in Vu, where £ is a 
phenomenological parameter. By the (anti) symmetry 
of reversible kinetic coefficients, this implies a piece 
—(piPjAfi in the stress cry, which we recognize as the 
active stress of Equation [8] 

Below is some recommended reading on recent appli- 
cations of active hydrodynamics t o the cytoskeleton. 
Active stresses have been shown jl!2| | to be impor- 
tant in the formation of the contractile ring around 
the equator of an animal cell about to divide. Oscilla- 
tions in experiments on non-adhering fibroblast cells 
have been understood through the interplay of active 
stresses in the cortical actin laye r with the entry of 
calcium through ion channels 113j . Strong depar tures 
from the fluctuation-dissipation theorem (FDT) |l!4j 
have been seen in a permanently crosslinked actin net- 
work studded with myosin motors, and analysed us- 
ing active hydrodynamics extend ed to the case of a 
permanently crosslinked network |l!5l Ill6l. Strongly 
non-FDT behaviour has been predicted 117] in an ac- 
tive hydrodynamic treatment of a stiff filament in an 
active medium, and conne ction s drawn to t he dynam - 
ics of auditory hair cells [1181 ] and axons 
The idea of actin-based active transport par allel to 
the cell membrane has been proposed in Ref. 121]. 
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swimming cells 
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FIG. 9: Relative viscosity of suspensions of chlamy- 
domonas (measured at shear rate 5 s _1 ) versus 
volume fraction cj>. Solid symbols: live cell data; 
crossed symbols: dead cell data. Reprinted fig- 
ure with permission from S. Rafai et al. 125] 
(http://prl.aps.org/abstract/PRL/vl04/i9/e098102)." 
Copyright 2010 by the American Physical Society. 



approach was moti vated by cytoskeletal physics, but 
more recently [1241 ] the focus has been on swimmers. 



G. Experiment on active suspensions 



1. Viscosity measurements on microbial suspensions 



2. From motors and filaments to active hydrodynamics 



Kruse and Jiilicher [122j . motivated by actomyosin 
phenomena, studied active stresses in one-dimensional 
microscopic models of motors and filament bundles. 
In an important series of pape rs, Liverpool, Marchetti 
and collaborators [H, [9ll Il23| have carried out a sys- 
tematic construction of the coarse-grained equations 
of motion of active filaments, for systems on a sub- 
strate as well as in bulk suspension, starting from a 
microscopic description of motors gripping filaments 
and moving along them. They obtain equations of 
the form propos ed b y p urely phenomenological the- 
ories [ll|, Ihil. Il6l. [l9ll42l l43j . but each parameter in 
the coarse-grained theories is traced to a microscopic 
process invol ving mot ors and filaments. An important 
outcome [HI, I9llll23| is phase diagrams including the 
various well-known instabilities of actin solutions such 
as bundling, with a connection to experimentally ac- 
cessible control parameters such as filament concen- 
tration and strength of motor activity. The initial 



The prediction [90j that activity alters the viscosity 
of suspensions, as discussed in sec tion IIV1 has re- 
cently been tested by Rafa'f et al. [l25l . Il26| on the 
motile alga chlamydomonas, a contractile swimmer, in 
a co nven tional rheometer, and by Sokolov and Aran- 
son [l27j ] on the bacterium Bacillus subtilis, an exten- 
sile swimmer, by monitoring the decay of an imposed 
vortex as well as by measuring the viscous torque on a 
rotating particle. The idea was to check whether the 
change in viscosity could be attributed directly to the 
swimming activity. Fig. [S] shows the viscosity of the 
microbe suspension as a function of concentration <j>, 
and compares the values obtained with live and dead 
cells, for <j) upto 25%. The dependence of viscosity on 
(j) was about twice as s trong for live as for dead cells. 
For B. subtilis in [127], the viscosity was suppressed 
by as much as a factor of 7 in some cases. These exper- 
iments not only confirm the predictions of Hatwalne 
et al. [9(| but also provide strong support for the idea 
that a medium suffused with active processes such as 
swimming should be thought of as a distinct type of 
material. 
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FIG. 10: Turbulence at zero Reynolds number in a sitting 
drop of Bacillus subtilis, viewed from below. Scale bar 
= 35 /jm. Reprinted figure with permission from C. 
Dombrowski et al, Phys. Rev. Lett. 93, 098103 (2004); 
http : //prola. aps . org/ abstract /PRL/v93/i9/e098103. 
Copyright (2004) by the American Physical Society. 



2. Order and instability of active suspensions 



There is ample qualitative evidence of the instability 
of active ordered suspensions predicted in reference 
and discussed in section IIVBI It would appear 
that collections of bact eria can't swim straight even 
if they set out to do so |l28l - ll30j |. and that their own 
active stresses are responsible. Bacterial turbulence is 
seen in Fig. [TUJ The role of instabilities in focusing the 
bacterial conce ntration h as been stressed by Goldstein 
and coworkers Il28l 1 13 ill . 



We have already seen one imitation of active matter 
in section IIII B 21 namely, a vertically agitated hori- 
zontal layer of bits of wire. Several efforts are under- 
way to make colloidal active matter, including beads 
or rods ha lf co vered with platinum and immersed in 
H 2 2 [H [l32j]. The platinum catalyzes 2H 2 2 -> 
2H 2 + 2 , and the asymmetry leads to directed mo- 
tion as a result of osmotic gradients. These hav e been 
analyzed theoretically in References jl03l Il33l Il34j . 
Elcctroosm qtic p ropulsion has been explored by Lam- 
mert et al. [135j A magnetically driven artificial flag - 
ellum has been demonstrated by Dreyfus et al. [136| . 
The study of self-propelled liquid drops on a surface, 
a common sight on a hot skillet in the kitchen, is a 
subfield in itself and cannot be reviewed in any de- 
tail here. The driving force is generally evaporation 
of one or more constituents from the liquid. A cou- 
ple o f recent examples are the work of Thakur et al. 
[137] on nematic drops and Chen et al. 138 1 on oil 
drops driven by surface tension gradients. It should 
also be possible to make motors by forcing asymmet - 
ric elastic dimers with nonequilibrium noise [l39l Il4dj | . 
An intriguing new theoretical development holds out 
the possibility of using mechanical agitation to make 
Stokcsian swimmers from collections of elastic dimers 
141]. In all these cases a well-known general principle 
[601 ] seems to work: polarity together with sustained 
dissipation leads to motility. 



3. The structure factor of fish shoals 



VI. CONCLUSION AND OUTLOOK 



A remarkable underwater acoustic measurement [107] 
of the static structure factor of the number density 
of fish shoals, for wavenumbers from 10 -4 to 10 -2 
m reveals a power law <~ fc~ 1,5 , and a rapid wave- 
like response of fish to perturbations, both of which 
are consistent with the Toner- Tu [ll| work described 
above. However, to deal with the fact that the fish are 
in water, a more complete theory of active suspensions 
at high Reynolds number, or at least in the unsteady- 
Stokes description, is clearly needed. Power count- 
ing on Equation 1101 shows that nonlinearities are rele- 
vant in dimension d < 4. An actual renormalization- 
group calculation for this problem is even more daunt- 
ing than that for the Toner- Tu model [ll[ and there- 
fore has not been done. A useful first step w ould be 
a comparison of the fish-shoal observations [l07| to 
the unsteady-Stokes predictions [l6T [ discussed under 
"Fish Schools" in section HVl 



What has the active-matter effort achieved? It has 
opened a new area of fundamental physics - the hy- 
drodynamics of liquid crystals and suspensions in an 
entirely novel regime. It is part of a larger process: 
biology is increasingly concerned with the mechanical 
forces exerted in living systems, and the active-matter 
approach provides the framework in which to ask these 
questions. It has made specific predictions, and exper- 
imenters are managing to check some of them. The 
impact on biology is already bei ng f elt, from the sub- 
cellular scale to starling flocks |l42| . In cell biology, 
tissue is clearly the next mo untain, a nd some steps on 
it have already been taken [l43l Il44l ] . The active hy- 
drodynamic framework, suitably adapted, should ap- 
ply to self-d riven par ticles in more complex media and 
geometries [I45l4l50j . 

The approach to active matter advocated in this re- 



view shares an important weakness with all coarse- 
grained theories: an abundance of phenomenologi- 
cal parameters. If thermal Brownian motion and 
active processes are present in comparable magni- 
tude, a clear distinction can probably be made be- 
tween parameters determined primarily by equilib- 
rium physics and arising strictly from activity. How- 
ever this nice separation is lost if noise an d co cfh- 
cients themselves are functions of activity [126j . as 
they probably are in bacteria and certainly are in 
granular matter. This limits somewhat the predic- 
tive power of coarse-grained theories. This is a limita- 
tion we're used to in equilibrium physics: Ginzburg- 
Landau free-energy functionals predict temperature 
dependences best near continuous phase transitions, 
not over an entire phase. In addition, living systems 
are finite in space and time and thus not ideally suited 
for gradient-expanded approaches. We must admit 
these limitations and accept humbly the need for mi- 
croscopic theories. 



There are two imperatives if the field is to make^a" 
real contribution to biology. Future theory work must 
aim to integrate the combined effects of growth, cross- 
linking, treadmilling, multiple species of filament, and 
the coupling of the cytoskeleton to the cell mem- 
brane [l2l| . Future experimental work must aim for 
quantitative tests of existing predictions of active- 
hydrodynamics, on controlled experimental systems. 
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